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1 Introduction 



The external field problem for matrix models has received recently much attention due 
to the work by Kontsevich |IJ. Generically, the (N x iV hermitian) one- matrix model in 
an external field is defined by the partition function 



Z[A;N] = JDX e Ntr{Ax-MX)) (L1) 



where Vq(X) is some potential. As follows from the results by Kontsevich, this model 
with a cubic potential Vq(X) oc X 3 describes Witten's formulation of 2D topological 
gravity. 

The partition function ( |1.1| ) can be calculated by the standard methods of solving 
matrix models. While the orthogonal polynomial technique can not be used due to the 
presence of the matrix A, the method || of Schwinger-Dyson equations written in terms 
of eigenvalues of A has been applied recently 0, ||| to explicitly solve the Kontsevich model 
by genus expansion. An analogous solution in genus zero has been obtained || for the 
case of the Kostov-Metha potential 

V (X) = ±X 2 -a\ogX, (1.2) 

when it is associated with an external field problem for the Penner model ||. 



The external field A in Eq.(l.l) plays a role of a source for correlators of trX fc which 
can be obtained by differentiating w.r.t. A and putting then A = 0. An alternative way 
to calculate these correlators is to introduce the most general potential 



OG 



V(X) = Y J tkX k (1.3) 

fc=0 

and to consider the partition function 

Z[t.;N) = [ DX e - try(x) (1.4) 



whose derivatives w.r.t. t^s reproduce the same correlators of trX fc 's provided one puts 
V(X) = NVq(X) after the differentiation. In this approach, the role of an external field 
is played by the set of couplings {tk}. It is evident that the two ways of introducing an 
external field are equivalent. As far as we know, this equivalence has not been elaborated, 
however, in the literature. 

In the present paper we prove the exact relation between the partition functions ( |1 . 1| ) 
and QDf): 

Z[A;N] = e% tT ^Z[t. ;aN], (1.5) 



which is valid provided A and {t } are related by the Miwa-type transformation 

t k = ltrA- k - f 8 k2 for n > 1, t = tr logA" 1 (1.6) 
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and N — > oo which makes t^'s to be independent variables. The proof is based on the 
fact that the Schwinger-Dyson equations for both models are equivalent. Based on this 
equivalence, we formulate then a method for calculating the partition functions order by 
order of genus expansion which is closed in spirit to that proposed by Gross and Newman 
Q for the unitary matrix model and for the hermitian one with a cubic potential. We 
calculate explicitly the partition function in genus one and apply the result to calculate 
genus one contribution to correlators of conformal operators with the puncture operator. 
Our genus one results support the conjecture that the model (|1.1| ) with the potential 
( |1.2|) (and therefore the model ( |1.4|) with the potential ( |1.3|) ) is associated with the c = 1 
case, if a is identified with the cosmological constant, in the same sense as the Kontsevich 
model is associated with c = 0. 



2 The Virasoro constraints 

A simplest way to prove the equivalence between the partition functions (|1.1|) and ( |1.4|) is 
to show that the Schwinger-Dyson equation for the model ( |1 . 1|) with the potential (1.2), 



° 2 +NA JL- aN A e~f trA2 Z[A] = 0, (2.1) 



[dA 2 8A 
is equivalent to the set of the Virasoro constraints, 

L n Z[t] =0 for n > -1 (2.2) 

where 

oo q n 8 8 

^ = £H fc ^ + £^-f-, (2.3) 

/%_0 ^^n+k k=0 k ^^n—k 

imposed on the partition function ( |1.4| ). 

Eq.( |2.1| ) has been discussed in Ref. [BJ. The extra gaussian factor in front of Z is 
introduced to cancel one arising from the integral (one can easily see this is the case of 
a — 0). The commutation of this gaussian factor with the derivatives makes the sign of 
the second term in parenthesises opposite to what would be naively expected. As for the 
Virasoro constraints (|2.2|) , (|2.3|) , they were obtained in Refs.Q using the invariance of the 
integral (|L|) under the shift X -> X + e n X n+1 . 

Now our purpose is to make a change of variables from (the eigenvalues of) A to {t} 
in Eq.( [2.1|) similarly to what was done in Refs.[|l(| |ll[] to derive the continuum Virasoro 
constraints from the Kontsevich matrix model (in that case only odd A;'s entered the 
Kontsevich-Miwa transformation ( |1.6| )). The only what we need is the chain rule 



8 00 8 



which is a consequence of Eq. (|1.6| ) . 
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Making this change of variables, one obtains from Eq.( |2.1| ) 



,\-n-2 r „-Ttr.\- 



A- n - 2 L n e-- trA Z = (2.5) 



with L n given by Eq. ( |2.3| ) provided that 

d 

—Z = -aNZ (2.6) 

which emerges formally for n = —2. The last equation coincides with the normalization 
condition for aN x aN hermitian one-matrix model. We have reproduces, therefore, 
exactly the Virasoro constraints (|2.2j ), providing the partition functions are related by 



Eq.Q. 

The fact that the Virasoro constraints (|2.2j ), fl2.3|) can be obtained from the Schwinger- 
Dyson equation for the Kontsevich model making the transformation (|1.6j ) with both even 
and odd /c's was noticed by Marshakov, Mironov and Morozov ||10|| . It was not recognized, 
however, that the proper external field problem is given by the partition function (|1 . 1|) 
with the potential ( |1.2| ). 



3 Loop equations 

Loop equations are successfully applied for studying the matrix models since the work by 
Kazakov |12| (for a review of recent results, see e.g. Ref.||13||). 

An insertion of the one-loop operator (with A being the Laplace transformed momen- 
tum which corresponds to a loop of the length I) is given by 

--E^- 1 ^. (3.1) 



This operator can be rewritten in terms of d/dA. Due to Eq.( |1.6| ), one gets 



Q 00 Q 

dXi k=0 1 dt k 



where Aj stands for an eigenvalue of A. As N — > oo when Eq.( p.2|) becomes exact, one 
can define the density of eigenvalues of A: 

1 N 

PW = YjH^-^) (3-3) 
iV i=i 

so that 

' 19 4 (3.4) 



^(A) NdX5p(X)' 
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Similarly, one can relate (the derivative of) the loop source V(A), which is defined by 
Eq.(|1.3|), to p. Using Eq.(|1.6|), one gets 



x 1 



V'{\) = £ kt k X k ~ 1 = £ t r - NX, (3.5) 

k=l j=l A J ~ A 



where the sum over k is convergent if A < min | A, |. As N — > oo Eq. (|3.5| ) determines V 
to be the Hilbert transform of p: 

This equation expresses unambiguously each term of the A-expansion of V (A) via p pro- 
vided the support of p{x) vanishes at some finite interval which includes the point x = 
(this is an analog of the above condition A < min | Aj |). All the integrals below are rigor- 
ously defined providing p(A) possesses this property. The discussion of general properties 
of Eq. (|3~6|) is beyond the scope of the present publication. 

We compare now two equations. The first one is the Schwinger-Dyson equation for 
the model ( (Lip written in terms of Aj pi 

92 +Ev 2 v^4)+"4-«» , )« f5i ' i %l = »' < 3 - 7 > 



<9Af f^i^i- Xj \9\ dX 3 



The second one is the loop equation written for the one-loop average 



which reads (see, e.g. JT3[] ) 

du V'(uj) 



IMX-ufM-WW + ^wm. (3.9) 

This equation is equivalent to the Virasoro constraints ( |2.2| ), ( |2.3| ) which can be obtained 
by expanding both sides of Eq. (|3.9| ) in 1/A. 

It follows immediately from Eqs.( |3.1|) , ( 3.4|) , (3-8|) that the r.h.s. of Eq.( [3.9|) coincides 
with the first term (with the double derivative) on the l.h.s. of Eq.( p.7|) . In order to 
compare the remaining terms, one inserts Eq. (|3.6|) into the l.h.s. of Eq. (|3.9| ) and calculates 
the contour integral over uj by taking residuals at the two poles. Applying again Eqs. ( |3.1| ) . 
(p.4|), (|3.8|), one gets exactly the remaining terms on the l.h.s. of Eq.(|3.7|) providing the 
normalization condition ( |2.6| ) is satisfied. This completes the proof of the equivalence 
between Eqs.(^.7|) and 



4 Comparison of genus zero solutions 

It is instructive to compare the known (one-cut) genus zero solution of loop equation (|3.9|) 
with that of Eq. (|3.7| ) which has been obtained in Ref.||. 
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The genus zero solution to Eq. ( |3.9| ), which was first obtained for a general V(A) that 
involves both even and odd powers of A by Migdal flH|, reads 



Wo (A) = / p ,^ + 4tA + 4c (4.1) 

ov ; J Cl Am (A - u) ^u 2 + Abu + 4c V ; 



with 6 and c given by 



dio V'(lj) f duj uoV'iu) 

0, / — -= V ; = 2aN. (4.2) 



Id 2-ni ^Juj 2 + Abu + 4c ' ic-j 27rz V^ 2 + 4ku + 4c 
Inserting Eq. (|3.6|) and doing the contour integral as before, one gets from Eq.( |4.1| 



\W Q {X) = \ 



N 



VX 2 + AbX + Ac-X-j 



dx 



p(x) N\ 2 + Ab\ + Ac 
x — X \ y/x 2 + Abx + 4c 



(4.3) 



while Eq.flO) yields 



\ I dx - P(X) + b = 0, \\dx , XP ^ + c-3b 2 = a + l. (4.4) 

2 ' Vx 2 + Abx + Ac 2 J Vx 2 + Abx + Ac 2 



One sees that ( [4.3| ), ( |4.4j ) coincide with the corresponding formulas of Ref. 0. The 
W(x) which was introduced there is related to our Wo by W(x) = jjWo(X) + \ while 
the definitions of b and c are the same. One more formula which would be useful for 
applications can be obtained by differentiating Eq.(fO|) (or (|4.1|) ) w.r.t. a: 

1 dW {X) _ 1 

N da VA 2 + 46A + 4c' 1 ' ' 

Being expanded in j, this expression reproduces Eq.(5.20) of Ref.||. An advantage of 
the approach of Ref . pi is that it allowed to calculate log Z itself. 

For an iterative solution of Eq. (|3.7|) which is considered below, we shall need an explicit 
expression for the irreducible two-loop correlator 



Applying Eq.( [3.4|) to Wo given by Eq.( |4.3|) , one gets in genus zero 
w(x v 1 Xp - 26(A + n) + 4c , , b 2 -c 

W ^,fl) 2(A _ ^2 V / A 2 + AbX + 4c ^2 + Abfi + 4c ' ^lA,Aj (A 2 + 46A + 4c) 2 

(4.7) 

which coincides with the known result by Ambj0rn et al. 0. 

The solution (|4.3|) , ( |4.4| ) is simplified if p(A) is a symmetric function p(A) = p{— A). The 
first equation in (|4.4|) yields then 6 = while only even powers of x enter the second one 
(as well as Eq.( [4.3| )). This case corresponds to the so-called reduced hermitian one-matrix 
model, i.e. to vanishing odd times t 2m +i- It implies for the Miwa transformation (|1.6|) 
that only even fc's are present. This situation is complementary to the original Kontsevich 



6 



model when only odd fc's appear. The loop equations for the reduced hermitian one-matrix 
model can not be formulated, however, entirely in terms of the even times. The odd times 
inevitably appear for this model to higher orders. This does not happen for the case of 
the complex matrix model (see, e.g. [JTBJ) which is formulated entirely via the even times. 
It is the model which is complementary in this sense to the Kontsevich model. A study 
of the external field problem for the complex matrix model will be published elsewhere. 



5 The genus one solution 

Our idea of how to solve Eq.( |3.7| ) (or ( |3.9| )) iteratively is similar to one which has been 
proposed by Gross and Newman Q for the unitary matrix model and for the hermitian 
one with a cubic potential. 

Let us introduce the new variables 

1 



B p = 2 p - 1 {2p-l)\\ J dXp(X) 



(A 2 + 46A + 4cf + ^ ' 
A 



C p = 2P-\2p-l)\\ fdXp(X) __( a+ i) 5p0 (5.1) 

J (A 2 + 46A + 4cf + 2 

where b and c are determined by Eq. fl4.4p which reads now 

B + b = 0, C + c - 3b 2 = (5.2) 

((— 1)!! = 1 by definition). Let us define F g — genus g contribution to logZ — by the 
formula 

oo 

log Z = J2N 2 - 29 F g . (5.3) 

9=0 

Now our conjecture is that F g would depend at 1 < g < oo only on B p and C p for p < P 
where P is some finite number (the larger g the larger P). This is in contrast to the 
t-dependence of F g which always depends on the whole set {£&}. Such a behavior of F g 



for the Kontsevich model has been advocated recently by Itzykson and Zuber [13 . 

Let us find explicitly the genus one correction to the genus zero results which are 
described in the previous section. Defining 

W l( X) = J^F, (5.4) 

with 5/5V(X) given by Eq.(|3.4j) and substituting into Eq.(p.7|), one gets the following 
linear equation for Wi(A): 



■^Wo(X) + A)Wi(A) + f dxp(x) Wl(X) x _ + ^ (A, A) = 0, (5.5) 



where W (X) and W (X, A) are given by Eqs.Q4.3p and ((4.7|), respectively. 
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For an arbitrary p, the solution of Eg. (fT5|) requires tedious calculations. Some sim- 
plifications occur if p is a symmetric function p(X) = p(— A). As is discussed above, this 
corresponds to the reduced hermitian matrix model (i.e. to an even potential V). For 
this case, one gets 6 = and B p = in Eqs.(|5.1|) and ( |5.2| ) which simplifies calculations. 
Using the property p(X) = p(— A), Eq.( |4.3|) can be rewritten as 



jfW (\) + A = VA 2 + 4c - j dxp(x) 



x 



VA 2 + 4c 



x 2 - A 2 Vx 2 + 4c' 



(5.6) 



By a direct differentiation of this formula one gets 



d_ 

dc 



■ N 



W (A) + A) 



2(1 - Gi) 
VA 2 + 4c 



(5.7) 



which will be extensively used below. We shall need as well the following rules of differ- 
entiation: 



dCp 
dc 
6d 
6V(X) 



cC 2 

i-c a 



Sc 



2c 



5V(X) (A 2 + 4c)l(Ci-l) : 



+ 



12c 



5 ) 
|2 



(5i 



(A 2 + 4c)f (A 2 + 4c)^ 

which are easy to derive from Eqs.( ^.l|) , (|5.2|) . 

Using the property p(A) = p{— A), let us represent the second term on the l.h.s. of 
Eq.( |5.5|) as the linear operator /C: 

AWi(A) -xWxix) 



KW X {\) = / dxp(x)K(X,x)W 1 (x) = fdxp(x 



A 2 — x 2 



(5.9) 



In matrix notations, when A and x are replaced by the eigenvalues Aj and Xj, this operator 
becomes aniV x JV matrix JC^i- Such an operator was considered in Ref.0. Eq.( p75|) can 
now be rewritten as 



jrW (X) + X)W 1 (X) + / dxp{x)K(X, x)W 1 (x) 



N (\2 



(A 2 + 4c) 2 



(5.10) 



It is the equation which is solved below. 



The form of the operator /C and the fact that Wi(X) = —Wi(—X) suggest the following 
ansatz 

oo A 

w 1 (x) = y: — n —^ (5.H) 



n+T 



n=0 (A 2 + Ac) 

where the coefficients {A n } are functions of {C p }. It is easy to calculate of how the 

operator fC acts on the 'basis vectors' 1/(A 2 + 4c) n+ 2. One first calculates for n = and 
then obtains a general formula by applying (d/dc) n . The result reads 



J dxp(x)K(X, 



■l) n 



x 



(* 2 + 4cH 2»(2n-l)! 



'cT 
dc 



2W (X) + X 



VX 2 + Ac 



(5.12) 



The term arising from the action of (d/dc) n on 1/a/A 2 + 4c equals to the first term on the 
l.h.s. of Eq.( |5.10| ) with the minus sign and cancels it when inserted into Eq. ( |5.10| ). The 
other terms can be easily calculated using Eq.( |5.7| ). Calculating the derivatives, one gets 
finally 

1 1 _ c n C i 

(2^ (A) + A + /C) L - r= M2 ■ ,\n ~E 



(A 2 + 4c) n + i (A 2 + 4c)n p f 2 2P-i(2p - 1)!! (A 2 + 4c)«- P +i ' 

(2Vy (A) + A + £) , } = = 1 for n = (5.13) 

V A 2 + 4c 

where the operator notation for the l.h.s. of Eq. (|5.10 ) has been used. This formula is 
similar to that by Gross and Newman j|] while the definition of the 'moments' C p is 
different for our model. 

One sees now that the r.h.s. of Eq.( |5.10| ) is reproduced by the n = 2 term so that 
Eq.( |5.10|) is satisfied if all A n = except 



C 2 C c 

Al = ^m-c 1 y> A2 = ^\~c[- (5 ' 14) 



Therefore, we have found the genus one solution to be 

^ l(A) = * {6(l-C 1 ) 2 (A 2 + 1 4c)l + l-V(A 2 + 4c)l} ■ (5 ' 15) 

It is worth noticing that we have obtained the explicit genus one solution to the reduced 
model with an arbitrary potential V(A) = V(— A). Previously it was calculated [|nj by 
solving loop equations only for the case of a polynomial potential with the highest power 
A 6 . We have verified that the result for the quartic potential can be reproduced by our 
formula ( |5.15| ) when p(X) is such that only £2 and £4 are nonvanishing so that 

Ci = 1 + 2t 2 - 24t 4 c, C 2 = 24t 4 , C p = for p > 3 (5.16) 

and Eq. ( |5.2|) for c reads 

- 2t 2 c + 12t 4 c 2 = a. (5.17) 

The general formulas which express the 'moments' C p via {t2k} for a polynomial V(A) 
look similar. 

Using the collection of formulas from Eq. (|5~8|) , one integrates Eq. ( p.!5| ) to obtain 

F 1 = ~]og{c{l-C 1 )}. (5.18) 



While our genus one result looks similar to that of Ref. 0, |T^], the coefficient in front of 
the logarithm is now ^ instead of h. As is discussed in the next section this is related 
to the fact that our matrix model is associated with c = 1 and not with c = 0. 
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6 The relation to c=l 



We discuss in this section of how our genus one results support the conjecture of Ref.0 
(which was based on the genus zero calculations) that the model ( |1 . 1| ) with the potential 
( |1.2| ) is associated with c = 1 providing a is identified with the cosmological constant. 
This model has been called the Kontsevich-Penner model. 

Let us start with the case of A — > oo when the Kontsevich-Penner model is reduced to 
the standard Penner model [§] which corresponds |J to b = 0, c = a, B p = and C p = 
for p > 1 in our notations. This case can be easily recovered by the genus one solution 
from the previous section. One gets from Eq. (|5.18|) 

F 1 = ~bga (6.1) 



which exactly reproduces the corresponding result of Distler and Vafa |L7j for the Penner 
model. Notice that we did not take the 'double scaling limit' to obtain this result since 
our model is associated with the continuum case. 

The correlators of conformal operators O n with the puncture operator V can be ob- 
tained for our model in genus one by applying d/da to Eq. fl5.15| ) and expanding in 1/A. 
The only C depends explicitly on a so that our collection of formulas ( |5.8| ) should be 
supplemented by 

dc 1 



da 1 — C\ 
By a direct differentiation of Eq.( |5.18D one gets 



(6.2) 



d 2 1 [1 C 2 C 3 2C| \ 

da* 1 ~ 12(1 - Ci) 2 [J + c(l - + (1 - Ci) + (1 - CO 2 J ( ' 

which corresponds to genus one contribution to the correlator of the two puncture op- 
erators V. One can combine this result with the genus zero calculation of Ref. [[]] to 
give 

" ; F = log(-C) (6.4) 



da 



where C is defined by 



C = c + * ^J l + 7 J^ + 77 ^ + 7 ^) (6.5) 



i2A"(i -df \c (i-co (i-co (i-cy 

The meaning of this formula is disscussed below. 

To calculate the next correlators, one differentiates Eq.( p,15| ) w.r.t. a. The result reads 

d f C 2 - cC 3 2cC 2 1 1 



N—Wi(X) 



da 1V ' i(l-Ci)3 (1-C!) 4 J 6(A 2 + 4c)l 

, f 1 2cC 2 ] 1 10c 1 , , 

Ul-CO 2 (1-^)3/(^ + 40)1 (l-C 1 ) 2 (A 2 + 4c)i' 1 ' ] 
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The explicit expressions for the correlators of 2 k with V, which can be obtained by 
expanding the r.h.s. of Eq. ( |6.6| ) in 1/A 2 , are rather involved while some simplifications 
occur after the shift (|6.5|) . As an example we present the sum of the genus zero and genus 
one results for the correlator of V and 2 (the dilaton operator in the case of the Penner 
model) which reads 



(0 2 V) 



o 2 



dt 2 da 



-2C 



1 



6A 2 (1 - Ci) 



2C 2 1 

:i - d) + c 



(6.7) 



When expressed in terms of {t 2 k}, this formula is to be compared with the corresponding 
correlator for c = 1 CFT. 



Let us speculate now on an nonperturbative essence of our results. Eq. Q6.4j) looks like 
the equation d 2 F/dt 2 ) = u which expresses the susceptibility u via F for the Kontsevich 
model while Eq. ( |6.5| ) is an analog of the string equation (to the given order of the genus 
expansion) which relates u to the 'times' {t^}- The transformation ( |1.6| ) looks like the 
transformation which has been introduced by Kharchev et al. |]TB| for the generalized 
Kontsevich models to relate their partition functions to the (reduced) KP r-function 
while differs from it by the fact that the role of the cosmological constant is played now 
not by t but by the 'conjugate' (in the sense of Eq.( [2.6|) ) variable a. As is shown in 
Ref.||18||, the whole description of c < 1 CFT's can be obtained in this framework and 



the set of Virasoro and W- constraints can be constructed. Our conjecture indicates, 



therefore, that the set of the Virasoro constraints (^2), (2J3) which is constructed with 
the use of the variable to — 'conjugate' to the cosmological constant a — can be used 
alternatively to describe a c = 1 CFT. 

One more argument in favour of our conjecture is a similarity between the Virasoro 
constrains ( |2.2| ), ( |2.3| ) and those for the case of the Witten's generalization [|19[ of the 
intersection theory on moduli space analytically continued to k — —3 (or to p = —1 in 
the notations of Ref.poj) when it is equivalent in a certain sense to the Penner model. It 



would be very interesting to study the relation between all these approaches. 
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